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Quantum Random Walks do not need a Coin Toss 
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Classical randomized algorithms use a coin toss instruction to explore different evolutionary branches 
of a problem. Quantum algorithms, on the other hand, can explore multiple evolutionary branches 
by mere superposition of states. Discrete quantum random walks, studied in the literature, have 
nonetheless used both superposition and a quantum coin toss instruction. This is not necessary, and 
a discrete quantum random walk without a quantum coin toss instruction is defined and analyzed 
here. Our construction eliminates quantum entanglement from the algorithm, and the results match 
those obtained with a quantum coin toss instruction. 

PACS numbers: PACS: 03.67.Lx 



I. MOTIVATION 

Random walks are a fundamental ingredient of non- 
deterministic algorithms , and are used to tackle a wide 
variety of problems — from graph structures to Monte 
Carlo samplings. Such algorithms have many branches, 
which are explored probabilistically, to estimate the cor- 
rect result. A classical computer can explore only one 
branch at a time, so typically the algorithm is executed 
several times, and the estimate of the final result is ex- 
tracted from the ensemble of individual executions by 
methods of probability theory. To ensure that different 
branches are explored in different executions, one needs 
non-deterministic instructions, and they are provided in 
the form of random numbers. A coin toss is the sim- 
plest example of a random number generator, and it can 
be included as an instruction for a probabilistic Turing 
machine. 

A quantum computer can explore multiple branches 
in a different manner, i.e. by using a superposition of 
states. The probabilistic result can then be arrived at 
by interference of amplitudes corresponding to different 
branches. Thus as long as the means to construct a vari- 
ety of superposed states exist, there is no a priori reason 
to include a coin toss as an instruction for a (probabilis- 
tic) quantum Turing machine. This is obvious enough, 
and indeed continuous time quantum random walks have 
been studied without recourse to a coin toss instruction 
0- Nevertheless, a coin toss instruction has been consid- 
ered necessary in construction of discrete time quantum 
random walks (see for instance, Refs.0,3). In this arti- 
cle, we demonstrate that this is a misconception arising 
out of unnecessarily restrictive assumptions. We explic- 
itly construct a quantum random walk on a line without 
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using a coin toss instruction, and analyze its properties. 

There also exists confusion in the literature about dif- 
ferent scaling behavior of discrete and continuous time 
quantum random walk algorithms (see again, Refs.0,Q), 
because the former have been constructed using a coin 
toss instruction while the latter do not contain a coin 
toss instruction. Our work eliminates this confusion in 
the sense that scaling behavior of discrete and continuous 
time quantum random walk algorithms, both constructed 
without a coin toss instruction, would coincide. There- 
after, a quantum coin would be an additional resource; if 
its inclusion can improve scaling behavior of some quan- 
tum algorithms, that should not be a surprise. 



II. QUANTUM RANDOM WALK ON A LINE 

A random walk is a diffusion process, which is gener- 
ated by the Laplacian operator in the continuum. To con- 
struct a discrete quantum walk, we must discretize this 
process using evolution operators that are both unitary 
and ultra-local (an ultra-local operator vanishes outside 
a finite range To begin with, consider the walk on a 
line. The allowed positions are labeled by integers, and 
the simplest translation invariant ultra-local discretiza- 
tion of the Laplacian operator is 

H\n) cx [_| n -l)+2|n)-|n + l)] . (1) 

The corresponding evolution operator is 

U(At) = exp(iHAt) = l + iHAt + 0((At) 2 ) . (2) 

With a finite At, U has an exponential tail and so it is not 
ultra-local. The evolution operator can be made ultra- 
local by truncation, say by dropping the 0((Ai) 2 ) part, 
but then it is not unitary. One may search for ultra-local 
translationally invariant unitary evolution operators us- 
ing the ansatz 
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but then the orthogonality constraints between different 
rows of the unitary matrix make two of {a, b, c} vanish, 
and one obtains a directed walk instead of a random walk. 

One way to bypass this problem and construct a ultra- 
local unitary random walk is to enlarge the Hilbert space 
and add a quantum coin toss instruction, e.g. 



U = 



£[|TXT 



This modification however brings its own set of caveats. 
If the coin state is measured at every time step (in other 
words, if the coin is classical), one gets no improvement 
over the classical random walk. With a unitary coin evo- 
lution operator, that entangles the coin state with the 
position state, the quantum walk performs better than 
the classical walk in certain algorithms. But in this case, 
the final results depend on the initial state of the coin, 
because quantum evolution is reversible and not Marko- 
vian. For example, the final state distribution of the 
quantum walk depends on whether the initial coin state 
was |t), or ||), or some linear combination thereof. To 
get around this initial coin state sensitivity, further algo- 
rithmic modifications such as averaging over initial coin 
states, or intermittent coin measurements, or use of mul- 
tiple coins, have been suggested, but they still leave a 
feeling of something to be desired. 



A. Getting Rid of the Coin 

The way out of the above conundrum is familiar to 
lattice field theorists @ . It has also been used to simulate 
quantum scattering with ultra-local operators , and to 
construct quantum cellular automata |8|]. In its simplest 
version, the Laplacian operator is decomposed into its 
even/odd parts, H = H e + H a , 
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The two parts, H e and H Q , are individually Hermitian. 
They are block-diagonal with a constant 2x2 matrix, and 
so they can be exponentiated while maintaining ultra- 
locality. The total evolution operator can therefore be 
easily truncated, without giving up either unitarity or 
ultra-locality, 



|® |n + l)(n| + |i)(||® |n-l)(n|l. (4) u ( At ) 



i(H c +H )At ifl t A( iH a At 



0((A*) 2 ) (8) 
U e (M)U (At) + 0((At) 2 ) . 



The quantum random walk can now be generated using 
U e U as the evolution operator for the amplitude distri- 
bution ip( n i t)i 



t/>(n,t) = [c/ e c/ ] 0) , 



(9) 



The fact that U e and U do not commute with each other 
is enough for the quantum random walk to explore all 
possible states. The price paid for the above manip- 
ulation is that the evolution operator is translationally 
invariant along the line in steps of 2, instead of 1. 

The 2x2 matrix appearing in H e and H D is propor- 
tional to (1 — ox ), and so its exponential will be of the 
form (cl -I- isci), |c| 2 + |s| 2 = 1. A random walk should 
have at least two non-zero entries in each row of the evo- 
lution operator. Even though our random walk treats 
even and odd sites differently by construction, we can 
obtain an unbiased random walk, by choosing the 2x2 
blocks of U e and U„ as Qj). The discrete quantum 
random walk then evolves the amplitude distribution ac- 
cording to 
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(11) 



U e U Q \n) = _[i| n -l) + | n )+i|n+l)-|ra+2(-l) n )]. (12) 



B. Relation to the Walk with a Coin 

Our construction of discrete quantum random walk has 
exchanged the up/down coin states for the even/odd site 
label. In the language of lattice field theory, this strat- 
egy resembles staggered fermions |fg, while that with a 
coin (or spin) is akin to Wilson fermions 0. Indeed, an 
explicit relation between our construction and that with 
a coin can be established. Let 

•<»■•>- us:; t) ) < i3 > 

describe the amplitude distribution in a two-component 
notation. Then Eas. H9ll2|) are equivalent to the evolution 
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Here, for clarity, we have denoted the basis states for 
iff by | TV). The symmetric coin operator C mixes the 
up/down components of iff. The walk operator U dis- 
tributes the amplitude equally between remaining at the 
same site and moving to the neighboring sites. The pro- 
jection operators (l±cr 3 )/2 pick out the amplitude com- 
ponents that move forward and backward. Finally, the 
operator o\ interchanges the up/down components of iff, 
producing what Ref. 10] has called the flip-flop walk. 

It is also instructive to note that while the diffusion 
operator H has the structure of a second derivative, its 
two parts, H e and H Q , have the structure of a first deriva- 
tive. This split is reminiscent of the "square-root" one 
takes to go from the Klein-Gordon operator to the Dirac 
operator. For a quantum random walk with a coin, this 
feature has been used to construct an efficient search al- 
gorithm on a spatial lattice of dimension greater than 
one [ToL Hl|. Reanalysis of that problem is in progress, 
without using a coin, in view of our results ^3] • 



III. ANALYSIS OF THE WALK 

It is straightforward to analyze the properties of the 
walk in Ea. i|12fl using the Fourier transform: 



^(M)^ e ' h ?i(M) 
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The evolution of amplitude distribution then follows 

i>{k,t) = e lWkt $+(k,0) + e- lWkt ^^(k,0) , (22) 

where ip±(k,0) are the projections of the initial ampli- 
tude distribution along e±. The amplitude distribution 
in the position space is given by the inverse Fourier trans- 
form of ijj(k,t). While we are unable to evaluate it ex- 
actly, many properties of the quantum random walk can 
be extracted numerically as well as by suitable approxi- 
mations. 

Consider a walk starting at the origin, ip (n, 0) = S n fi. 
Its amplitude distribution at later times is specified by 
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27r J -ir \/l + sin 2 k 

— i sincjfct sin k + cos u>kty/l + sin 2 k 
i sin ujf-t 



This walk is asymmetric because our definitions treat 
even and odd sites differently We can get rid of the 
asymmetry by initializing the walk as ij) s (n, 0) = (<5 n) o + 
S n .i)/ v2. The walk is then symmetric under n <-> (1— n), 
and the amplitude distribution evolves according to 



^s,±(fc,0) = 



±1 



2 J 2(1 + sin 2 ft) 



(25) 



in ft ± yl + sin 2 ft 
1 + e lk sin k ± e lk \/l + sin 2 k 



The evolution of the amplitude distribution in Fourier 
space is easily obtained by splitting it into its even/odd 
parts: 



r [ J ) . ra-.l) [.W,-) '<•(/.-.(!) 



M(k) 



ie tk smk i cos ft 
i cos ft ie~ lk sin ft 



(18) 



(19) 



The unitary matrix M has the eigenvalues, A± = e ±luJk 
(this ± sign label continues in all the results below), 



2tt 



dk e 



- ikn 



2(1 + sin 2 ft) 



is\TLUJkt(e lk — sin ft) + cos Wktyl + sin 2 ft 



(26) 



y i sinujkt(l + e lk sin ft) + coswfct e %k \J 1 + sin 2 ft J 

Figs. 1-2 show the numerically evaluated probability dis- 
tributions, after 32 time steps, for asymmetric and sym- 
metric quantum random walks respectively. Note that, 
by construction, the distributions after t time steps re- 
main within the interval [— 2t + 1, 2t]. 



X± = sin 2 ft ± i cos ft V 1 + sin 2 ft , uik = cos ^sin 2 ft) , 

(20) 

with the (unnormalized) eigenvectors, 



e± ex 



sin ft ± y/l + sin 2 ft 



sin ft ± vT + sin 2 ft 



(21) 



A. Asymptotic Behavior of the Walk 

For large t, a good approximation to the distribution, 
Eq. (|2"6*|) , can be obtained by the stationary phase method, 
as in Ref. p| . The integral is periodic, and a sum of 
terms of the form 



I(n 
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FIG. 1: Probability distribution after 32 time steps for the 
asymmetric quantum random walk ip . 




FIG. 2: Probability distribution after 32 time steps for the 
symmetric quantum random walk ip s . The dark curve denotes 
the smoothed distribution of Eo. l|39ft , 



The highly oscillatory part of the integrand is determined 
by 4>(k, n,t) — — kn± bJ^t, while the remaining part g(k) 
is bounded. Simple algebra yields 



~dk 

d 2 L0 k 
dk 2 

d 3 uj k 
dk 3 



2 sin k 



\J 1 + sin 2 k 
2 cos k 
~ (1 + sin 2 fc)3/2 ' 
4sinfc(l + cos 2 k) 



(1 



fc) 5 / 2 



The stationary point of the integral, k = 

n 2 sin fco 

a = — = T 

1 VI + sin 2 fc 



(28) 
(29) 
(30) 

has to satisfy 
(31) 



which has a solution only for n S [— V^t, V2t]. 
We now separately consider the three cases: 

(1) \n\ > \[2t: There is no stationary point in this case. 
For \n\ = (v2 + e)i, \d(f>(k)/dk\ > e, and repeated inte- 
gration by parts shows that the integral falls off faster 
than any positive integer power of et. 

(2) |n| = y/2t: In this case, there is a stationary point 



of order 2 at fco = =Fsgn(n)7r/2. The integral is therefore 
proportional to i -1 / 3 . Explicitly, 



ip s (V2t,t) = ct- 1 ' 3 

ip s (-V2t,t) = ct- 1/3 
1 



( 1 + ^) cos (tI) 




, R r - w 0.355 
27T3 1 / 6 V3. 



(3) \n\ < y/2t: There are two stationary points in this 
case, fcoi G (— 7r/2,7r/2) and fco2 = 7r — fcoi, with 



sin fcg = T 



V4P 



~2 ' 
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V4t 2 - 2n 2 (At 2 - n 2 ) 



(35) 



(36) 



The integral is therefore proportional to t 1 / 2 . In terms 
of the phase, 



0o = —koin + uj ko t - (7r/4) 
the distribution amplitude is 
1 

1p8 = 



(37) 
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V4t 2 - 2n 2 /(2t + n) 



(38) 



VAt 2 -2n 2 /^At 2 -n 2 
((1 - i)n + 2t) / (2t + n) 



The smoothed probability distribution, obtained by re- 
placing the highly oscillatory terms by their mean values, 
is 



I "0s I smooth 



At 1 



TrVAt 2 - 2n 2 [At 2 - n 2 ) 



(39) 



(Here, the n <-> (1 — n) symmetry can be restored by 
replacing n by (n— 5).) As shown in Fig. 2, it represents 
the average behavior of the distribution very well. Its 
low order moments are easily calculated to be, 



V2t 



-y/2t 



|'0 



s I smooth^ 71 



n=-V2t 
s/2t 



\n\ ■ \ip. 



s I smooth 



dn — t , 



(40) 
(41) 



n 2 \A\Loot b dn = 2(2-V2)t 2 . (42) 

n=— V2t 

The following properties of the quantum random walk 
are easily deduced from all the above results: 
• The probability distribution is double-peaked with 
maxima approximately at ±v2t. The distribution falls 
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off steeply beyond the peaks, while it is rather flat in 
the region between the peaks. With increasing t, the 
peaks become more pronounced, because the height of 
the peaks decreases slower than that for the flat region. 

• The size of the tail of the amplitude distribution is 
limited by (ei) -1 ~ < -1 / 3 , which gives An> = A(ei) = 
0{t l / z ). On the inner side, the width of the peaks is 
governed by \ J' k t\-' l l' 2 ~ i" 1 / 3 . For |n| = (-y/2 - 5)t, this 
gives An < — A(St) = Olt 1 ^). The peaks therefore make 
a negligible contribution to the probability distribution, 

• Rapid oscillations contribute to the probability distri- 
bution (and hence to its moments) only at subleading 
order. They can be safely ignored in an asymptotic anal- 
ysis, retaining only the smooth part of the probability 
distribution. 

• The quantum random walk spreads linearly in time, 
with a speed smaller by a factor of \/2 compared to a 
directed walk. This speed is a measure of its mixing 
behavior and hitting probability. The probability distri- 
bution is qualitatively similar to a uniform distribution 
over the interval [— \/2t, \/2t\. In particular, m th moment 
of the probability distribution is proportional to t m . 

These properties agree with those obtained in Ref. 0] 
for a quantum random walk with a coin-toss instruction, 
demonstrating that the coin offers no advantage in this 
particular set up. (Extra factors of 2 appear in our re- 
sults, because a single step of our walk is a product of 
two nearest neighbor operators, U e and U a .) It is im- 
portant to note that the properties of the quantum ran- 
dom walk are in sharp contrast to those of the classical 
random walk. The classical random walk produces a bi- 
nomial probability distribution, which in the symmetric 
case has a single peak centered at the origin and variance 
equal to t. 



B. The Walk in Presence of an Absorbing Wall 

The escape probability of the quantum random walk 
can be calculated by introducing an absorbing wall, say 
between n = and n = — 1. Mathematically, the absorb- 
ing wall can be represented by a projection operator for 
n > 0. The unabsorbed part of the walk is given by 



i/j(n,t + l) = P n >oU e U V>(?M) 

1 



(43) 



U e U ip(n,t)- - J n) _i(*V>(0,*)-V(M)) 



with the absorption probability, 



n>0 



(44) 



Fig. 3 shows the numerically evaluated probability distri- 
bution, in presence of this absorbing wall, after 32 time 
steps and with the symmetric initial state. Compari- 
son with Fig. 2 shows that the absorbing wall disturbs 



the evolution of the walk only marginally. The probabil- 
ity distribution in the region close to n = is depleted 
as anticipated, while it is a bit of a surprise that the 
peak height near n — increases slightly. As a result, 
the escape speed from the wall is little higher than the 
spreading speed without the wall. As shown in Fig. 4, we 
find that the first two time steps dominate absorption, 
^ s ,abs(i - 1) - 0.25 and P s , abs (i = 2) = 0.375, with very 
little absorption later on. Asymptotically, the net ab- 
sorption probability approaches P s ,abs(oo) ~ 0.4098 for 
the symmetric walk. (We also find, for the asymmetric 
walk starting at the origin, P ,abs(oo) « 0.2732.) This 
value is smaller than the corresponding result P a b s (oo) = 
2/ir for the symmetric quantum random walk with a coin- 
toss instruction 

Thus the part of quantum random walk going away 
from the absorbing wall just takes off at a constant speed, 
hardly ever returning to the starting point. Again, this 
behavior is in a sharp contrast to that of the classical 
random walk. A classical random walk always returns to 
the starting point, sooner or later, and so its absorption 
probability approaches unity as t — ► oo. 
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FIG. 3: Probability distribution after 32 time steps for the 
symmetric quantum random walk, with an absorbing wall on 
the left side of n — 0. 
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FIG. 4: Time dependence of the probability for the symmetric 
quantum random walk to remain unabsorbed, in presence of 
an absorbing wall to the left of n — 0. The walk gets within a 
few per cent of the asymptotic escape probability in just two 
time steps. 
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C. Comparison to the Walk with a Coin 

The above results bring out the differences of our 
quantum random walk construction compared to that of 

Refs.[llD3: 

(1) We have absorbed the two states of the coin in to 
the even/odd site label at no extra cost. This is possible 
because, due to its discrete symmetry, the walk with a 
coin effectively uses only half the sites (e.g. for a walk 
starting at origin, the amplitude distribution is restricted 
to only odd sites at odd t and only even sites at even t). 
Our walk makes use of all the sites at every instance. 

(2) It can be seen from Eas. i|15l4|) that, at every time 
step, ^ has 50% probability to stay put at the same lo- 
cation, while the walk with a coin has no probability to 
remain at the same location. Yet, both achieve the same 
spread of amplitude distribution, as exemplified by the 
moments in Eas. H4QI42[) . This means that our walk is 
smoother — more directed and less zigzag. 

(3) When the coin is considered a separate degree of 
freedom, quantum evolution entangles the coin and the 
walk position. On the other hand, when the coin states 
are made part of the position space, as we have done, 
entanglement disappears completely — only superposition 
representing the amplitude distribution survives |15| . 
This elimination of quantum entanglement would be a 
tremendous advantage in any practical implementation 
of the quantum random walk, because quantum entangle- 
ment is highly fragile against environmental disturbances 
while mere superposition is much more stable. The cost 
for gaining this advantage is the loss of short distance 
homogeneity — translational invariance holds in steps of 
2 instead of 1. 



IV. EXTENSIONS AND OUTLOOK 

The quantum random walk on a line is easily converted 
to that on a circle by imposing periodic boundary condi- 



tions. When the circle has N points, the only change in 
the analysis is to replace the integral over k in the inverse 
Fourier transform by a discrete sum over fc-values sepa- 
rated by 2tt/N. Since the quantum random walk spreads 
essentially uniformly, there is not much change in its be- 
havior. All that one has to bear in mind is that, on a 
long time scale, unitary evolution makes the walk cycle 
through phases of spreading out and contracting towards 
the initial state. 

Going beyond one dimension, the strategy of con- 
structing discrete ultra-local unitary evolution operators 
by splitting the Hamiltonian in to block-diagonal parts is 
applicable to random walks on any finite-color graph • 
One just constructs 2x2 block unitary matrices for each 
color of the graph, and the single time step evolution 
operator becomes the product of all the block unitary 
matrices. In particular, the d— dimensional hypercubic 
walk can be constructed as a tensor product of d one- 
dimensional walks, using 2d block unitary matrices. 

Our results clearly demonstrate that discrete quantum 
random walks with useful properties can be constructed 
without a coin toss instruction. The addition of a coin 
toss instruction may still be beneficial in specific quan- 
tum problems. A coin is an extra resource, and there 
are known instances where the addition of a coin toss in- 
struction makes classical randomized algorithms have a 
better scaling behavior compared to their deterministic 
counterparts 0. One may hope for a similar situation 
in the quantum case too, keeping in mind that a careful 
initialization of the quantum coin state would be a must 
in such cases. 

A clear advantage of quantum random walks is their 
linear spread in time, compared to square-root spread in 
time for classical random walks. So they are expected to 
be useful in problems requiring fast hitting times. Some 
examples of this nature have been explored in graph the- 
oretical and sampling problems (see Refs.0, Q for re- 
views), and more applications need to be investigated. 
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